In [1]:
# Standard setup block for running Python code
from __future__ import division, print_function
import os
if os.path.split(os.getcwd())[-1] == "Lab notebooks":
os.chdir("../../")
print("Moved to experiment root directory")
from Modules.processing import *
from Modules.plotting import *
plt.style.use("Config/plotstyle.mplstyle")
%matplotlib inline
Got to the lab about 12:30 PM.
Homed turbine, tow, y-, and z-axes.
Tank depths is 2.44 m -- good.
Did dummy tow at 1 m/s to mix up tank.
12:51 PM -- Added a cup of seeding and did another dummy tow.
Did third dummy tow and Vectrino measurements look good.
1:02 PM -- Starting Wake-1.0-0.375. Reset 9188 chassis first.
Vectrino is at y=-1.4781 m, z=0.3026 m. -- Good.
4:07 PM -- Added a cup of seeding before run 30.
5:23 PM -- Finished Wake-1.0-0.375. RPOS(z) at bottom limit is 0.03600 m. RPOS(y) at negative limit is -1.5248 m.
5:29 PM -- Starting Wake-1.0-0.5.
Noticed that the Vectrino data from Wake-1.0-0.375 was a smaller file than the rest. Checked it out and the Vectrino stopped recording early, cutting off data that is not important. The tow itself lines up nicely with the linear encoder measurements.
6:09 PM -- Stopping Wake-1.0-0.5 for now.
In [91]:
Section("Wake-1.0-0.5").process(nproc=1, nruns="new")
In [92]:
WakeProfile(1.0, 0.5).plot_mean_u()
watermark()
In [89]:
r = Run("Wake-1.0-0.375", 3)
plt.plot(r.time_vec_all, ts.smooth(r.u_all, 30))
plt.plot(r.time_ni_all, r.tow_speed_all)
plt.show()
In [90]:
# Setting things up
setup_time_sec = 3600
n_setup_times = 4
# Seconds per tow at each speed
seconds = {0.8 : 90 + 200,
1.0 : 84 + 240,
1.2 : 80 + 300}
# Wake runs
runs_per_wake = 43
nwake = 2.8
wake_speed = 1.0
# Runs for strut torque and strut covers
strut_covers_runs = 22
strut_torque_runs = strut_covers_runs*2
# Tare drag and torque
tare_drag_runs = 14
tare_torque_runs = 30
tare_drag_dur = 200
tare_torque_dur = 120
total_time_sec = n_setup_times*setup_time_sec + runs_per_wake*nwake*seconds[wake_speed] \
+ strut_covers_runs*seconds[wake_speed] + strut_torque_runs*tare_torque_dur \
+ tare_drag_runs*tare_drag_dur + tare_torque_runs*tare_torque_dur
total_hours = total_time_sec/3600
hours_per_day = 10
total_days = total_hours/hours_per_day
print(total_hours, total_days)